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ABSTRACT 



The mass function of galaxy clusters is an important cosmological probe. Differences in the selection method could potentially lead 
to biases when determining the mass function. From the optical and X-ray data of the XMM-Newton Follow-Up Survey, we obtained 
a sample of galaxy cluster candidates using weak gravitational lensing, the optical Postman matched filter method, and a search for 
extended X-ray sources. We developed our weak-lensing search criteria by testing the performance of the aperture mass statistic on 
realistic ray- tracing simulations matching our survey parameters and by comparing two filter functions. We find that the dominant 
noise source for our survey is shape noise at almost all significance levels and that spurious cluster detections due to projections of 
large-scale structures are negligible, except possibly for highly significantly detected peaks. Our full cluster catalog has 155 cluster 
candidates, 116 found with the Postman matched filter, 59 extended X-ray sources, and 31 shear selected potential clusters. Most of 
these cluster candidates were not previously known. The present catalog will be a solid foundation for studying possible selection 
effects in either method. 

Key words. Gravitational lensing - Galaxies: clusters: general - X-rays: galaxies: clusters - Catalogs - Surveys 



1. Introduction 

Because the dynamical or evolutionary timescale of clusters of 
galaxies is not much shorter than the Hubble time, they retain a 
'memory' of the initial conditions for structure formation (e.g., 
iBorgani & Guzzol200lb . The population of clusters evolves with 
redsh ift, and this evo lution depends on the cosmological model 
(e.g.. lEke et al.|[T996h : therefore, the redshift dependence of the 
cluster abundance has been used as a cosmological test (e.g., 
IVikhhnin et al.ll2003l;lHenrvll2004 . The dependence of the clus- 
ter abundance on cosmologi cal parameters can be o btained ei- 
ther from analytical models jPress & Schechtei 1974.) or, mor e 
reliably, from A^-body simulations (e.g., iJenkins et al.l 1200 Ih . 
What such models do predict is the abundance of dark matter 
halos as a function of redshift and mass. 

Clusters can be selected by various m ethods: optical, X- 
ray emission, weak lensing (e.g.. ISchneider, 1996.) . and - using 
future surveys - the Sunyaev-Zeldovich effect (SZE). Optical, 
X-ray, and SZE selection of clusters depends on the bary- 
onic content of clusters, which - compared to the predicted 
dark matter density - is a minor fraction of the clusters' con- 
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stituents. Optical selection depends on the star formation history, 
and X-ray detection require s a hot intracluster medium (ICM). 
IWeinberg & Kamionkowskil (|2002|) predict that up to 20% of all 
weak lensing clusters have not heated their ICM to a level de- 
tectable with current X-ray missions. Searching for clusters with 
SZE is a very promising method but has yet to produce first sam- 
ples. 

Common to all methods except X-ray selection is that they 
are prone to projections along the line-of-sight. Spectroscopy is 
an essential tool in distinguishing real clusters from chance pro- 
jections on the sky. X-ray emission, which depends on the square 
of the local gas density, is not easily affected by line-of-sight pro- 
jections but is susceptible to other sorts of biases, e.g., heating 
of the ICM in mergers. 

All four selection methods are sensitive in different redshift 
regimes. Optical and X-ray searches depend on the luminosity 
distance of clusters. SZE is nearly redshift independent. Weak 
lensing surveys typically cover the redshift range of 0.15-0.7. 
Clearly, no cluster selection method is ideal, and understanding 
their biases and limitations is important for precision cosmology 
using clusters. 

Several galaxy clusters have already been found using weak 
gravitational lensing in recent years: in the FORSl co smic 
shear survey (iMaoli et al.l 120011; iHetterscheidt et al.ll2005l) . one 
mass peak clearly coincides with an overdensity of galaxies. 
Dietrich et al. (2005) found a clus ter in the back ground of the 
super-cluster system A 222/223. IWittman et al.l (|2006), have 
published 6 new clusters detected with weak lensing from their 
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Deep Lens Survey. These examples clearly demonstrate that this 
method of cluster detection in fact works. 

The selection of clusters of galaxies by weak lensing, how- 
ever, also faces significant methodological challenges. Even be- 
fore the first spectroscopica lly confirmed w e ak-lensing detected 
cluster was reported ( Wittman et al.ll200lh . lErben et al.1 (|2000|) 
reported a highly significant tangential a fignment of galaxies 
arouii d an empty spot on the sky (see also lvon der Linden et al.l 
120061) . Two m ore of these dark clumps have been reported in 
the literature (Umetsu & Futamase 2000; Miralles et al. 2002!, 
but see also E rben et al. 20 03). Recently S chirmer et al. (20071) 
have published a catalog of 158 shear selected peaks, 72 of them 
associated with bright galaxy concentrations. Of course, the re- 
ality of these dark clump detections has to be considered with 
caution, given that even one of them would have a profound im- 
pact on our understanding of the evolution of dark and baryonic 
matter in the Universe. 

Several theoretical studies have recently shed some light 
on the problem of dark clump detectio n in w eak lensing 
surveys. Among th e m are iHamana et al. I ( 12004 H04) and 
iHennawi & Spergell (l2005l HS05) who have both used ray- 
tracing simulations through A^-body simulations to study the ef- 
ficiency, also called purity by other authors, and completeness 
of the detection of clusters of galaxies in weak lensing surveys. 
An important result of HS05 is that the efficiency, even in the 
limiting case of no intrinsic galaxy ellipticity, does not increase 
beyond 85%. The remaining 15% of shear selected peaks are due 
to projections of the large-scale structure along the line of sight 
and will be seen as dark clumps. These could very well account 
for the dark clumps studied so far in detail in the literature. The 
efficiency naturally drops further if more realistic noise caused 
by the ellipticity of the background galaxi es is assumed. The 
completeness was studied in more depth bv lH04i who find that, 
even with low significance thresholds in the selection process 
shear, selected samples will be incomplete, except at the highest 
masses. 

In this paper we describe a search for gala xy clusters in the 
public XMM-Newton Follow-Up Survey (Dietrich et al.ll2006l) 
and our private extension to the survey. We combine the results 
of three independent selection methods using the aperture mass 
statistic (Schneider 199^ for weak lensing sele ction, the optical 
matched filter algorithm (Postm an et al.|[l996l hereinafter P96) 
for optical selection of galaxy clusters, and a search for extended 
X-ray emission in the XMM-Newton data on our survey fields. 
The combination allows us to dig deeper into the mass function 
than we could do with weak lensing selection alone. 

This paper is organized as follows: Section |2] briefly de- 
scribes the XMM-Newton Follow-Up Survey, the optical data, 
and their reduction. Section |3] gives an overview of the X- 
ray data, their reduction, and then describes how we detected 
galaxy cluster candidates as extended X-ray sources. In Sect. |4] 
we present our implementation of the matched filter technique 
of |P96l and the resulting catalog of cluster candidates. We 
briefly summarize the aperture mass methods ( jSchneideiill996l) 
in Section |5] and develop our selection criteria used for cluster 
detection later in that section based on realistic ray-tracing sim- 
ulations. We discuss and summarize our findings in Sect. |6] The 
catalogs of cluster candidates are available in electronic format 
at the CDS. 

Throughout this work we assume an Q.^ - 0.3, Oa = 0.7, 
Hq = 70/i7okms~' cosmology and u se standard lensing notation 
(e.g.. lBartelmann & Schneideiil2001l) . 



2. The XMM-Newton Follow-Up Survey 

The XMM-Newton Follow-Up Survey (XFS) consists of two 
parts, a public and a private one. Both were conducted with the 
Wide-field Imager (WFI) at the ESO/MPG-2.2m telescope on 
La Silla, Chile. The survey provides optical imaging on fields 
for which deep, public XMM-Newton exposures exist. 

2.1. Public and private survey 

The public survey (ESO Program Id. 170.A-0789, PI. J. Krautter 
as chairman of the ESO working group for public surveys) was 
carried out in the framework of the ESO Imaging Survey (EIS) 
as a collaboration between ESO, the XMM-Newton Survey 
Science Centre (SSC), and a group at the Institut fiir Astrophysik 
und Extraterrestrische Forschung (lAEF) at the University of 
Bonn. The aim of the public survey, whose observations were 
concluded in November 2005, was to provide optical counter- 
parts to serendipitously detected X-ray sources in the southern 
hemisphere. To provide data for a minimum spectral discrimi- 
nation and photometric redshift estimates, the public survey was 
carried out in the B-, V-, R-, and I-passbands down to a limiting 
magnitude of 25 mag in the AB system in all bands. All data of 
the public survey taken before October 16, 2003 were reduced, 
calibrated, and publically released in July 2005 dDietrich et al.l 
2006). The total public survey comprises 15 WFI fields. Out of 
these, 4 are galactic fields and thus unsuitable for galaxy clus- 
ter searches. One field is a mispointing without X-ray data and 
high galactic absorption. The remaining 10 fields cover about 
~2.8sq.deg. in BVRI. 

The private extension of the survey (Program Ids. 70. A- 
0529, 71.A-0110, 072.A-0061, 073.A-0050, PI. P Schneider) 
with the goal of a weak lensing search for galaxy clusters has 
been conducted as a collaboration between the lAEF and the 
Astrophysikalisches Institut Potsdam (AIP). Originally targeted 
to observe 17 additional fields in B- and R-band and to be fi- 
nalized by October 2003, observations were concluded only in 
September 2006 due to weather and scheduling problems. Out 
of the 14 additional fields observed in R-band, our primary band 
for the lensing analysis, 13 are used for this work. This corre- 
sponds to all data obtained until September 30, 2005. B-band 
observations, proposed to allow for a red-sequence cluster search 
(iGladders & Yeell2000h . are available for 9 of fliese fields. 

A cluster survey carried out on these fields is obviously not 
an unbiased survey since it includes 5 XMM-Newton observa- 
tions that targeted known galaxy clusters. The impact on the to- 
tal number of clusters detected, however, is small, at least for the 
X-ray and optically detected clusters. The influence on the weak 
lensing sample is greater because weak lensing requires clusters 
to be fairly massive to be detectable with 2 m class telescopes. 
Previously known clusters are likely to be relatively massive and 
influence the smaller sample of weak lensing detected clusters 
more strongly. 

2.2. Optical data reduction 

We reduced the optica l XFS data using the publically avail- 
able GaBoDS pipeline (lErben et al.ll2005h with the Guide Star 
Cat alog ver s ion 2. 2 (GSC-2.2) as the astrometric reference and 
the IStetsonI (l2000h catalog for photometric calibration. A sub- 
set of our reduce d data was compared against the publically re- 
leased XFS data dDietrich et al. 2006) and found to be in very 
good agreement with this independent reduction. 
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We would like to point out that we applied the fringing re- 
moval procedure to all our R-band images to remove the low 
level fringing pr esent in WFI R -band. This is a difference to the 
data released bv iDietrich et al.l ilOOQ . Weight images describ- 
ing the relative noise properties of each pixel were created to 
mask bad pixels or columns, cosmic rays, and other image de- 
fects masked manually, like satellite tracks and ghost images due 
to internal reflections. 

Table [T]lists all coadded R-band images with center coordi- 
nates of the survey fields used in the present work, their 5cr lim- 
iting magnitude computed in the Vega system in a 2" diameter 
aperture, their seeing, as well as the effective (unmasked) area 
used for catalog creation (see Sect. 14. It , and the number den- 
sities of the resulting lensing (see Sect. I5.2l i and matched filter 
(Sect.l4!2]i catalogs. 

Two independent observations of the field RBS 0864 are 
used in the XFS. The WFI observations in the V- and R-bands 
done by Schindler et al. (RBS 0864-N in Table [B are centered 
on coordinates slightly north west of the galaxy cluster and were 
taken with a seeing constraint of V.'2. The re-observation in 
the B- and R-bands two years later in the course of the XFS 
(RBS 0864-S) is centered on the cluster and was observed with 
the seeing constraint of I'.'O of this survey. However, most data 
was taken well outside the specified constraint. Consequently, 
we used only the RBS 0864-N pointing for our weak lensing 
analysis. We will describe the full XFS data set including BVl 
passband data in a forthcoming paper (Dietrich et al., in prepa- 
ration). 

3. X-ray detection 

The X-ray luminosity of the hot intracluster gas scales with the 
square of its density. Thus, X-ray detections of clusters of galax- 
ies are relatively insensitive to projection effects. The high sen- 
sitivity of XMM-Newton allows us to find clusters out to very 
high redshifts. Two of the three most distant and spectroscopi- 
cally confirmed clusters were found sere ndipitously with X MM- 
Newton. These are at redshift z = 1.39 (i Mums et allllob S) and 
z = 1.45 ([Stanford et al.ll2006h . the latter located in the public 
XFS field LBQS 2212-1759. 

3.1. X-ray data reduction 

The archival XMM-Newton data of the XFS were reduced with 
the latest version of the Science Analysis System available at that 
time (SAS-6.5.0). We briefly describe the standard processing 
steps employed to create X-ray source lists from the observation 
data files. Appendix lAl gives an overview of all X-ray observa- 
tions considered in the present work. Some fields were observed 
more than once. Their data reduction requires additional steps 
described at the end of this section. 

SAS was used to generate calibrated event lists, exposure 
maps describing the spatial variations of the instruments' sensi- 
tivity, and images in five energy bands ranging from 0. 1-0.5 ke V, 
0.5-l.OkeV, 1.0-2.0keV, 2.0-4.5keV, and 4.5-12.0keV. We 
used the standard SAS flags #XMMEA_EM and #XMMEA_EP to filter 
the event lists when generating the images. For the FN detector 
we imposed additional restrictions on the allowed charge pat- 
terns. For the two lowest energy bands we accepted only single 
events; for the three highest energy bands only events not de- 
positing charge in more than two neighboring pixels were kept. 
We also excluded the energy range of 7.8-8.2 keV from the PN 
data to avoid the complex of Ni-Ka, Cu-Ko', and Zn-Ka fluores- 
cence lines of the detector and surrounding structure. 



Using the 15 images generated with evselect (5 images 
per camera and 3 cameras), the mask images, and the vignetted 
exposure maps, eboxdetect in local mode was run to generate 
a first source list. The minimum likelihood for a detection in 
this step was set to 5. The sources were excised from the image 
by esplinemap to model the background of the images. The 
task esplinemap allows to describe the background either by a 
two-component model (vignetted astrophysical and unvignetted 
particle background) based on ray-tracing of the instruments or 
by a 2-d spline with a user-defined number of nodes. For the 
XFS the decision about which approach to use was based on 
visual inspection of the data. 

With these background maps, eboxdetect was run in map 
mode to create another source list. The minimum likelihood for 
source detection was set to 4 in this step. This list was given 
together with the science images, masks, and background maps 
to emldetect, which performs a simultaneous maximum like- 
lihood multi-source PSF fitting in all energy bands. The free 
parameters emldetect fits are the source position, source ex- 
tent (core radius of a j6-model), and source count rate in each 
energy band and instrument. Derived parameters are the total 
source count rate, total likelihood of detection and likelihood of 
detection per energy band, likelihood of source extent, and four 
hardness ratios between the input energy bands. In our reduction 
we let emldetect fit up to two sources to one source position 
reported by eboxdetect. The minimum likelihood for a detec- 
tion was set to 6, and the minimum extent likelihood for a source 
to be considered as an extended source to 4. 

It should be pointed out that the result of the PSF fitting per- 
formed by emldetect cannot be better than the available PSF 
models. Especially at large off-axis distances the EPIC PSF is 
not very well known. 

For the two deepest combined fields (LBQS 2212-1759, 
250 ks, PG 1 1 15+080, 220 ks) the source detection process was 
slightly modified with respect to the other observations. In deep 
observations small inaccuracies of the background fit can lead to 
many spurious detections of extended sources. Since in the soft- 
est band (0.2-0.5 keV) both the MOS and PN cameras show spa- 
tial variations of the detector background and since the hardest 
band (4.5-12 keV) is strongly dominated by background, only 
the 3 bands in the range 0.5^.5 keV were used for source detec- 
tion in the two very deep observations. 

The source positions in the emldetect catalog have statis- 
tical errors on the order of l"-2", plus a systematic error due 
to an uncertainty in the attitude of the spacecraft, which has the 
same size. The latter can be corrected by cross-correlating the 
X-ray source positions with the more accurate positions of opti- 
cal sources. The task eposcorr was used to carry out this cross- 
correlation with the R-band catalogs on all fields. The production 
of the optical catalogs is described in detail in Sect. 14. 11 

Areas of the EPIC field-of-view that are dominated by bright 
extended emission of the XMM-Newton target were excluded 
from the survey by masking out a circular area around the target 
before starting the detection process. 

Multiple observations of one field are in general taken at dif- 
ferent roll angles and - since the center of rotation is not the 
center of the FOV - different positions. The first step in combin- 
ing observations is thus to bring them all to a common nominal 
position. This was done with the SAS task attcalc, which com- 
putes sky coordinates for event files. From these new event lists 
FITS images, exposure maps, and masks were created for the 
individual observations as described above. Science images and 
exposure maps were coadded weighted by the masks of the re- 
spective observation. The masks themselves are combined with 
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Table 1. Effective area and number density of galaxies of WFI fields in the XFS used in the cluster search. 



Field a 6 Seeing Area Number Density (WL) Number Density (MF) 

(J2000.0) (J2000.0) (mag) (arcsec) (arcmin-) (arcmin-^) (arcmin--) 



BPM 16274 


00: 


:50: 


:03, 


,2 


-52: 


:08: 


:17 


25, 


,38 


QSO BO 130-403 


01: 


:33: 


:01, 


,9 


-40: 


:06: 


:28 


24, 


,94 


CFRS 3h 


03: 


:02: 


39. 


,2 


+00: 


:07: 


:31 


25, 


,46 


RXJ0505. 3-2849 


05: 


:05: 


:20, 


,0 


-28: 


:49: 


:05 


25, 


,58 


RBS 0864-N 


10: 


:21: 


:03, 


,8 


+04: 


:26: 


:23 


25, 


,53 


QSOB1033-033 


10: 


:36: 


:23, 


,7 


-03: 


:43: 


:20 


25, 


,59 


SDSS J104433.04-012502.2 


10: 


:44: 


:33. 


,0 


-01: 


:25: 


:02 


25, 


,44 


MS 1054.4-0321 


10: 


:56: 


60, 


,0 


-03: 


:37: 


:27 


25, 


,35 


HE 1104-1805 


11: 


:06: 


:33, 


,0 


-18: 


:21: 


:24 


25, 


,47 


PG 1115+080 


11: 


:18: 


:17, 


,0 


+07: 


:45: 


:59 


25, 


,48 


CD -33 07795 


11: 


:29: 


:27, 


,2 


-34: 


:19: 


:55 


25, 


,04 


T Leo 


11: 


:38: 


:27, 


,1 


+03: 


:22: 


:10 


25, 


,47 


IRAS 12112+0305 


12: 


:13: 


:46, 


,1 


+02: 


:48: 


:41 


25, 


,51 


LBQS 1228+1116 


12: 


:30: 


:54, 


,1 


+ 11: 


:00: 


:11 


25, 


,09 


NGC 4666 


12: 


:45: 


08, 


,9 


-00: 


:27: 


:38 


25, 


,43 


QSO B 1246-057 


12: 


:49: 


:13, 


,9 


-05: 


:59: 


:19 


25, 


,21 


Field 864-1 


13: 


:41: 


:22, 


,4 


+00: 


:23: 


:52 


25, 


,45 


Field 864-9 


13: 


:44: 


:36, 


,0 


-00: 


:24: 


:00 


25, 


,56 


A 1882 


14: 


:14: 


:39, 


,9 


-00: 


:19: 


:07 


25, 


,73 


MKW9 


15: 


:32: 


:29, 


,3 


+04: 


:40: 


:54 


25, 


,56 


LBQS 2212-1759 


22: 


:15: 


:31, 


,7 


-17: 


:44: 


:05 


25, 


,52 


NGC 7252 


22: 


:20: 


:44, 


,8 


-24: 


:40: 


:42 


25, 


,57 


PHL 5200 


22: 


:28: 


:30, 


,4 


-05: 


:18: 


:55 


24, 


,92 



Total/average 



logical or. The source extraction continues on the combined im- 
ages, exposures maps, and masks as described above. 



3.2. Catalog of extended X-ray sources 

Catalogs of cluster candidates were generated from the 
emldetect source list. Sources with a detection of likelihood 
> 15 and extent likelihood > 4 were considered as potential 
cluster candidates. Any remaining large-scale inhomogeneities 
in the background are sometimes detected as spurious sources 
and the best-fit model of emldetect is often one whose ex- 
tent reaches the maximum value of 20 pixels (80"), hence only 
sources with an extent < 20 pixels were kept in the final cata- 
log. The extended X-ray sources found in this way were visu- 
ally screened and any obvious artifacts, often due to out-of-time 
(OOT) events or remaining background structure, were manually 
rejected. The optical images were visually inspected for possible 
counterparts of extended sources. In this step extended sources 
obviously associated with nearby galaxies were removed from 
the catalog. Grades were assigned to the quality of a cluster de- 
tection based on visual inspection of the optical and X-ray im- 
ages. Cluster candidates with grade "+" are obvious real clus- 
ters, often the ones one would select by eye. Extended X-ray 
sources graded with "o" are possible clusters but not as promi- 
nent as those graded with "+". This grade is assigned to clusters 
without a very obvious optical counterpart but a reliable X-ray 
detection. In some cases these may be systems at very high red- 
shift that are just barely visible in the optical images. Clusters 
graded "-" appear to be unreliable in the optical and X-ray im- 
ages, but were not rejected as obvious spurious sources based on 
their visual impression in the X-ray images. 

The full X-ray cluster catalog is available in electronic for- 
mat from the CDS. Cluster candidates in this catalog follow the 
naming scheme BLOX JHHMM.m+DDMM.m, where BLOX is 
the lAU registered acronym for "Bonn Lensing, Optical, X-ray" 



,88 


1011 


15.6 


21.6 


,01 


925 


7.8 


12.7 


,85 


965 


12.1 


17.0 


,83 


981 


18.2 


25.9 


,01 


948 


13.0 


20.2 


,78 


1037 


16.8 


22.0 


,73 


997 


16.9 


22.6 


,69 


1022 


17.7 


23.5 


,83 


1030 


13.4 


22.4 


,92 


1051 


12.3 


19.6 


,91 


927 


6.1 


19.1 


,74 


1005 


18.6 


23.9 


,88 


1028 


13.9 


21.9 


,01 


1014 


8.4 


14.3 


,87 


1074 


10.7 


17.8 


,80 


869 


13.0 


19.6 


,88 


1058 


16.8 


23.9 


,84 


1045 


16.9 


25.6 


,72 


1019 


18.6 


29.5 


,91 


1032 


10.7 


19.6 


,99 


970 


12.6 


19.9 


,70 


997 


20.3 


26.4 


,06 


1068 


7.8 


11.7 




23073 


14.1 


20.9 



detected cluster candidates. In the following we give comments 
on individual fields when appropriate. 



- RX J0505. 3-2849 - Two extended X-ray peaks are each 
found on both previously known RX J clusters in this field. 
In the case of RX J0505. 3-2849, we only list one; the other 
is probably a confusion with a double point source; 

- RBS 0864 - Two reductions of this field were done. One 
reduction excluded the target cluster from the analysis, while 
the other one was made with the cluster remaining in the 
data. A large number of spurious detections associated with 
OOT events were manually rejected. 

- MS 1054.4-0321 - Two extended X-ray sources were de- 
tected on the target cluster at about the same distance from 
the optical center of the cluster. We list both X-ray sources 
in the catalog. 

- HE 1 104-1805 - The only extended X-ray source found in 
this field coincides with a bright star. 

- LBQS 1228+1 116 - The X-ray data on this field is strongly 
affected by background flares. No clusters were found in the 
remaining shallow data. 

- MKW 9 - The target cluster dominates the center and several 
extended sources are detected in the cluster region, some of 
them on bright and large galaxies. 

- QSO B 1246-057 - The calibration problem that prevented 
the X-ray dat a on this field f rom being included in the public 
data release dPietrich et al.ll2b 06) were solved and the field 
could be included in the cluster search. 

- A 1882 - The X-ray image of this cluster shows three very 
extended so urces c entered around the nominal position of 
the cluster of lAbelll ( 1958). These regions of extended X-ray 
emission were excluded from the analysis. 

- NGC 7252 - No extended targets were found in this field. 
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4. Optical matched filter detection 

Clusters of galaxies can be optically selected in a multitude of 
ways, either from one passband alone or by combining color in- 
formation from two or more passbands with the positional in- 
formation on galaxies used in all methods. A review of a large 
number of optical detection methods has recently been published 
bylGal (2006). 

For this work we cho se the matched filter detection algorithm 
of iPostman et al.l (Il996l) . The |P96i method was selected because 
it is well-tested and efficient, it works on single passband cata- 
logs, and can thus be used for the entire area of the XFS. More 
elaborate detection schemes using multi-color information will 
be employed on XFS data in subsequent work (Dietrich et al. in 
preparation). 

4.1. Optical catalog creation 

The starting point for any optical cluster-search method is a cat- 
alog of galaxies. It is thus only prudent to discuss the creation of 
such catalogs in some detail before turning to a short description 
of the matched filter algorithm and how it was implemented for 
this survey. 

Unfortunately, WFI images are especially prone to internal 
reflections, producing prominent reflection rings around all sat- 
urated stars with blooming in the core. The brightest stars show 
more than one reflection ring, with increasing sizes and offsets 
from the position of the source. Additionally, extended halos and 
diffraction spikes occur around bright sources. The catalog cre- 
ation tries to mask regions aff'ected by these problems and do so 
as automatically as possible, but it still requires a large amount 
of manual masking. 

The catalog production uses SExtractor jBertin & ArnoutsI 

11996 ^ and starts with a very low signal-to-noise ratio (SNR) cat- 
alog. The weight images produced by the GaBoDS pipeline are 
used in all steps of the catalog creation. The sole purpose of this 
first catalog is to identify regions that should be masked. Masked 
regions will be passed on to SExtractor in a FLAG image. 

All objects of the initial catalog are put into cells of a grid 
whose size is chosen such that the average number of objects 
per cell is 2. This grid is smoothed with a Gaussian kernel 
with an FWHM of 1.4 pixels. We call the smoothed array S. 
Automatically adjusting the grid size such that a fixed number 
of galaxies per cell is reached on average allows one to keep 
the size of the Gaussian, which is necessary to reach the desired 
SNR in the object density distribution, as fixed in the program. 
The main advantage of this approach is that the values of the 
kernel can be stored in a matrix (fixed to a 5 x 5 array in the pro- 
gram). The convolution of the density grid with the kernel matrix 
is computationally much faster than re-evaluating the Gaussian 
kernel at each grid cell. In the next step the dynamic range S is 
limited. Every pixel with a value > 1.5 is set to 2. Every pixel 
with a value < 1 is set to zero; these pixels will be masked be- 
cause they are either in very low SNR parts of the image, such as 
the edges, or they are covered by extended objects, such as large 
foreground galaxies or very bright stars. The resulting array is 
called D. 

Additionally, any rapid change of object density in D, such 
as seen at the edges of reflection rings, is detected with a Sobel 
edge detection, i.e., an array containing the absolute values of 
the gradient of the array D. The gradient computation is imple- 
mented as a convolution of D with the two 3x3 convolution 
kernels that correspond to finite second-order, two-sided differ- 



entiation. Every pixel exceeding a threshold in D is flagged in 
the output image. 

Finally, the FLAG array, containing only values of 1 for pix- 
els to be flagged and for all other pixels, is smoothed again 
with the 1 .4-pixel FWHM Gaussian to account for the fact that 
the initial smoothing shrinks the areas not covered by objects. 
Every pixel with a non-zero flag value will be flagged in the out- 
put FLAG image, which is expanded to the size of the original 
WFI image from which the catalog was created. 

FLAG images created through this procedure reliably mask 
extended objects, bright stars, and the most prominent reflection 
rings. Exceptionally empty regions on the sky are only rarely 
masked erroneously. However, fainter reflection rings and stars 
of intermediate magnitude must still be masked by hand. Files 
describing the regions masked manually - either by circles or 
polygons for more complex shapes - must be generated by the 
user with tools such as the image viewer DS9. From them and 
the automatically generated FLAG image, the final FLAG image 
used in further steps in the catalog production is created. 

The image seeing is determined from a high-SNR catalog. A 
histogram of the FWHM of aU objects with a'3 < FWHM < 3'.'0 
is created and the image seeing is set equal to the FWHM of 
the bin with the most objects. The seeing is used as input for 
SExtractor's star-galaxy classifier CLASS_STAR. Since we are 
only interested in galaxies and not in stars, reliably separating 
them in the science-grade catalogs is important for galaxy cluster 
searches. 

Finally, two science-grade catalogs are created with 
SExtractor. Their only difference is the filter with which 
the detection image is convolved. The catalog for the optical 
matched filter search is made using a Gaussian kernel with an 
FWHM of 4 pixels. This kernel ensures that relatively few spu- 
rious detections of faint objects are made, but it has a lower com- 
pleteness at the low SNR end with a tendency to miss very small 
objects. Another catalog with SExtractor's default filter - a 
3x3 pixel pyramidal kernel - is created as the starting point for 
the weak-lensing catalog creation (see Sect.|5]l. 

4.2. The Postman matched filter catalog 

The Postman matched filter algorithm is described in detail in 
P96. The main features of this algorithm separating it from other 
single-band cluster detection schemes are (1) it uses photomet- 
ric and not only positional information (2) the contrast of struc- 
tures that approximate the filter shape with the background is 
maximized, (3) redshift and richness estimates of cluster candi- 
dates are produced as a byproduct, (4) the algorithm performs 
well over a wide range of redshifts. The main disadvantage is 
that a particular radial profile and luminosity function must be 
assumed. Clusters deviating from the expected profile will be 
detected only at lower significances or s upp ressed. Our imple- 
mentation follows the description of lP96l and l01sen et alj(ll999l 
099 hereinafter) with some modifications of the selection pa- 
rameters as outlined below. 

A series of programs and shell scripts is used to go from the 
object catalogs described in Sect. |4.1| to a catalog of cluster can- 
didates. The input catalog is filtered such that only objects with 
a high probability of being galaxies are kept. To this end only 
objects with a CLASS_STAR value < 0.5 are kept. At the faint end 
the catalog is cut at a magnitude lOcr above the sky background 
as measured in a 2" aperture to ensure a high completeness of 
the catalog. Objects brighter than 17 mag are also filtered. These 
are often bright stars, which are not correctly identified as such 
by SExtractor and would cause a serious contamination of the 
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input catalog. If any of these objects with R < 17 are bright, 
nearby cluster members, the matched filter signal of the cluster 
is decreased. However, the fainter cluster galaxies usually still 
lead to a significant detection of the galaxy cluster but with a 
redshift bias. This bias is introduced because omitting the bright 
cluster galaxies modifies the luminosity function. This redshift 
bias introduced by the cut on the bright galaxies is marginal, at 
least for the redshift range considered here. 

Likelihood maps S (x, y) are computed on a grid and saved as 
FITS images for a series of redshifts. We compute S with a spac- 
ing in redshift of Az = 0.1, starting at z = 0.1 and up to the last 
output grid in which the m* of the fiducial cluster model shifted 
to that redshift is not fainter than the limiting magnitude of the 
input catalog. For the XFS R-band observations, this is typically 
up to z = 0.9. Following the prescription of 099, we assume 
a fiducial cluster model at z = 0.6 with a Schechter luminosity 
function with = -21.63 mag, faint end slope a - -1.1, and 
a Hubble profile with physical core radius 7?c =0.1 /ivoMpc in- 
tegrated out to the cutoff radius rca = lO^c- We keep the pixel 
size of the output grid constant at Q.Spix/R^ at z = 0.6 but vary 
the area required for a detection in the output images with red- 
shift. The ^-correction needed for the redshift dependence of m* 
is computed for an elliptical galaxy with no evolution. 

Having created the likelihood maps, we are now faced with 
the problem of identifying clusters of galaxies in them; i.e., 
we have to decide which peaks in the likelihood maps are re- 
liable candidates for galaxy clusters. Peaks are detected with 
SExtractor. The main SExtractor detection parameters were 
adopted from 099 and are set as follows: 

- The minimum area for a detection MIN_AREA scales with the 
redshift and corresponds to nr^, where is the size of core 
radius projected on the sky; 

- The detection threshold DETECT.THRESH is set to 2; 

- No deblending of peaks is performed (DEBLENDJIINCONT^ 
1), so that all contiguous pixels above the detection threshold 
are counted as one cluster candidate. 

- A global background is estimated from the image. 

Catalogs from the individual output grids at different red- 
shifts are then matched by position using the LDAC program 
associate. Peaks present in at least 3 output grids with a min- 
imum significance of at least 3.5cr in one of them are kept as 
reliable cluster candidates. These parameters are slightly differ- 
ent from the ones adopted by 099, who used a threshold of 3cr, 
a detection in at least 4 output grids, and a minimum value of the 
richness parameter A^. The reason for different selection criteria 
are the different redshift regions of interest in the study of l099l 
and the work presented here. While 099 are chiefly interested in 
high redshift clusters, our search mainly aims for intermediate 
redshift clusters accessible with weak lensing, i.e., most of these 
clusters will b e at re dshifts z - 0.2 . . .0.3. Several of the crite- 
ria adopted bv l099l are biased against the search for clusters at 
intermediate redshifts. 

The significance of a cluster signal drops sharply once we 
look at output grids at redshifts higher than the cluster's redshift 
because the power law cutoff of the luminosity filter strongly 
suppresses signals from lower redshifts. For clusters at lower 
redshifts, the number of output grids at redshifts lower or equal 
to the cluster redshift is small. The requirement of 099 that a 
cluster candidate be detected in at least 4 output grids is thus 
heavily biased against the detection of clusters at the redshift 
interval we are especially interested in. We therefore relaxe the 
requirement on the number of output grids in which a cluster 
must be detected to 3, and have the lowest redshift output grid at 



z = 0.1 instead of z = 0.2. To compensate for the higher num- 
ber of spurious detections caused by this less stringent cut, the 
minimum significance to be reached in at least one output grid 
is increased from 3 to 3.5. 

We do not make any cuts on the richness parameter Ad. 
Clusters at higher redshifts need to be much more luminous than 
clusters at lower redshifts to be detectable. Consequently, the 
average richness parameter of clusters at higher redshift is much 
larger than that of lower redshift clusters. Imposing a minimum 
value for Ad would create a strong bias against low and interme- 
diate redshift systems. 

The selection criteria adopted by 099'have proven to be very 
successful by a high rate of spectroscopic ally confirmed clus- 
ters found in the EIS wide survey (e.g., [Hansen et al.1 120021; 
lOlsen et al ] |2003HBenoist et al.ll2002l) and we are thus confident 
that our slightly modified criteria are also successful. However, 
we have to note that the simulations on which these criteria were 
develope d are not able to predict the number of spurious de- 
tections. l099l attempted to simulate a pure background popu- 
lation of galaxies by randomizing the position of galaxies while 
keeping their magnitude fixed. This randomized population has 
a much smoother distribution than the field population of galax- 
ies. The reason is that the distribution of field galaxies is not a 
random field with white noise, but instead shows clustering and 
correlation. 

More recently, lOlsen et al.1 (l2007l O07) have used a 
more realistic approach by creating backgrounds that are cor- 
rela ted within bins of 1 mag following the prescription of 
Son eira & Peebles ! (Il978h . Their selection parameters are a min- 
imum significance of 3.5cr, minimum area of nr^, and number 
of redshift slices > 1, which are very similar to ours. 

We applied the matched filter algorithm to R-band catalogs 
from 23 fields of the XFS. The resulting catalog is available 
in electronic format from the CDS. It lists 116 candidate clus- 
ters. Thirteen of the systems in this catalog have been previously 
found by other authors and have either spectroscopic or photo- 
metric redshift information. All cluster candidates were visually 
inspected and graded as for to the grades given to the X-ray de- 
tected clusters. Cluster candidates with grade "H-" are obvious 
real clusters, often the ones one would select by eye, or they 
have and extended X-ray source as a counterpart. Matched-filter 
peaks graded "o" are probable detections but less prominent than 
those graded "-h" and may be more prone to projection effects. 
Candidates with a "-" grade are most likely artifacts or maybe 
in some cases very poor clusters or groups. 

Out of the total of 1 16 matched-filter peaks 49, were graded 
"+", 42 "o", and 25 received a grade of There is a clear 
correlation between the significance of a cluster candidate and 
its grade. The average SNR of a "+" rated cluster is 5.5cr and in- 
clude the highest SNR detections; "o" graded cluster candidates 
have an average significance of 4.9, while the unlikely candi- 
dates rated with "-" have an average value of (r,nax = 4.4. 

In the field RBS 0864, two clusters were detected 
independently on both pointings on this field. These 
are BLOX J1023.6-)-041 1.1 (RBS 0864 itself) and 
BLOX J1022.9-H0411.9. Their properties reported in the 
table of matched-filter detected clusters at the CDS are averages 
of both independent detections. 

A comparison of the matched-filter estimated redshifts of the 
13 previously found systems with their redshift values found in 
the literature is provided by Fig. \T\ Seven of these 13 clusters 
have spectroscopic redshift information, while the remaining 6 
cluster have only photometric redshift estimates. We point out 
that the redshifts of the photometric sample are not photometric 
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Fig. 1. Comparison of redshift estimates obtained from the 
|P96l matched filter with literature values. Solid circles denote 
data points with spectroscopic redshift information; open trian- 
gles are photometric redshift estimates from the literature. The 
matched filter underestimates the redshift as is shown by the 
best-fit lines to the spectroscopic sample (solid line) and the fuU 
sample (dashed line). 



redshifts in the classical sense of, e.g., template fitting but are the 
redshift estimators of other cluster finding methods. This is most 
n otably the redshif t estimator of the "cut-and-enhance method" 
of iGoto et al.l (12002 ). which provides 4 redshifts in this sample. 
The mean offset from the zero line for the spectroscopic sample 
is marginally significant with (Az/(1 + Ztme)) - -0.06 + 0.05, 
while the mean offset of the whole sample (Az/(1 + Zu-ue)) - 
-0.06 ± 0.07 is consistent with zero, where the error is simply 
the standard deviation of the scatter. The relatively small devia- 
tion from zero hides a significant bias of the matched filter red- 
shifts towards lower redshifts. The solid line in Fig. [T] is the best 
linear fit to the spectroscopic data points. The line is described 
by Az/(1 + Ztrue) = (0.00 ± 0.03) + (-0.16 ± 0.06)ztrae. Thus, the 
deviation from the ideal relation is significant at the 2.7cr level. 
The bias increases marginally if the photometric data points are 
included in the analysis. The best-fit line is then described by 

AZ/(1 + Ztme) = (0.07 ± 0.05) + (-0.29 + 0.12)Zt„e. 

Since the clusters were found on different fields with inde- 
pendent photometric calibration, errors in the photometric cali- 
bration can be excluded as the source of this systematic differ- 
ence. Possible sources of the underestimated redshifts are (1) the 
^-correction, which depends on the adopted galaxy model and 
(2) the luminosity function of the fiducial cluster model, specif- 
ically the value of M*. The discrepancy increases with higher 
redshifts. Our flux filter was constructed under the assumption 
of no evolution and thus stronger discrepancies are indee d ex- 
pected at higher redshifts. This has already been noted by lP96l 
who also found the estimated redshifts of higher redshift clusters 
to be systematically too low by about Az = 0.1 . . .0.2. 

The O07 paper carried out a matched-filter cluster 
search in the Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS) Deep. The CFHTLS field D4 coincides with 



Fig. 2. Histogram of the estimated redshift distribution of the 
matched-filter detected clusters. The open histogram represents 
all cluster candidates, the gray histogram the "+" and "o" rated 
candidates, and the black histogram only the "-h" rated clusters. 



our field LBQS2212-1759. We briefly compare the re- 
sults of these two independent cluster searches. They 
find 6 matched-filter selected cluster candidates inside the 
WFI FOV. Two of their detections, CFHTLS-CL-J221500- 
175028 and CFHTLS-CL-J221537-174533, coincide with 
matched-filter cluster candidates, BLOX J22 15. 0-1 750. 5 and 
BLOX J2215.6-1745.5, detected in this survey. Both clus- 
ter candidates received the best grade in either survey. 
The matched filter redshift estimates agree for CFHTLS-CL- 
J221500-175028/BLOX J2215.0-1750.5, while the 'O07 put 
CFHTLS-CL-J221537-174533/BLOX J2215.6-1745.5 slightly 
higher at z = 0.4 th an our estimate of z = 0.3. Two other cluster 
candidates oflOOTi coincide with possible clusters we found with 
other methods. CFHTLS-CL-J221620-173224, an "A" rated po- 
tential cluster at an estimated redshift of z = 0.7 matches our 
weak lensing detection (see Sect. |5]l BLOX J2216. 3-1733.0. 
CFHTLS-CL-J221606-175132, a "B" rated cluster candidate 
with a redshift estimate of z = 0.4 matches our weak lensing 
detection BLOX J22 16. 1-1751. 7. Our matched-filter cluster c an- 
didate BLOX J2214.4-1728.1 at z = 0.2 is not found bvlO07l 



5. Weak lensing detection 

The tidal gravitational field of a galaxy cluster causes a coher- 
ent tangential of the sheared images of background galaxies. A 
quantitative measure for this alignment, the aperture mass Map, 
was developed in Schneider (1996). There, Map is defined as a 
weighted integral over the dimensionless surface mass density k. 



-I ' 

•JsupU 



d^0 u{§)K{e) 



(1) 



where U{d-) - U{\9 - 9o\) is a radially symmetric weight func- 
tion with zero total weight to avoid the mass-sheet degeneracy 
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dSchneider & Seitj|1995h . The aperture mass in Eq. ([TJ is de- 
fined in terms of the surface mass density, but it is possible to 
find an expression that allows computing Map in terms of the 
observable shear y: 

Map(eo)= r d20e(,?)yt(0;0o), (2) 
where we define the tangential shear relative to a point 0o by 

r,(e;0o) = -1R[r(^ + 0o)e-'"^] . (3) 

Here, ()?, ^) are polar coordinates with respect to and the 
weight function is related to by 



(4) 



On real data the aperture mass can be estimated by a sum over 
the A^g galaxy ellipticities inside the aperture. 



1 

n ^ 



(5) 



Here n is the number density of faint background galaxies 
(FBG), and the tangential component of the ellipticity has been 
defined in analogy to Eq. ([Sjl. 

The SNR of a peak in maps of aperture mass can be esti- 
mated by using the fact that (Map) = holds in the case of no 
lensing. Then the RMS dispersion of the Map-statistic becomes 



= ^<Map>2, 



which is 



1/2 



usmg 



(6) 



(7) 



What remains to be fixed is the shape of the weight func- 
tion. If only white noise caused by the random ellipticity of 
background galaxies is present, i.e., if the noise is described by 
Eq. (|6l), the filter function Q should follow the shear profile of 
the cluster as closely as possible to increase the SNR of a galaxy 
cluster detection. The aperture mass then becomes a matched fil- 
ter technique for weak lensing detections of galaxy clusters. 

A^-body simulations predict the shape of collapsed halos to 
follow an NEW pr ofile. The shear of an NEW profile can be com - 
puted analytically jBartelmannll 1 9961: IWright & Brainerdll2000l) . 
but the resul ting expressions are c omplex and time-consuming 
to evaluate. Schirmer et al.l (l2007l) have proposed a simple ap- 
proximation to the NEW shear profile. 



2nfw(x) - 



tanh 



(8) 



that is much faster to compute than the full expression. The fore- 
going equation is not used directly for the computation of Map, 
but exponential cut-offs are multiplied to Eq. (O as — > and 



1 tanh(jc/xc) 



(9) 



The purpose of these cut-offs is (1) to avoid finite field effects as 
the weight function ([8j formally extends to infinity, but data is 
only available on a finite field; and (2) to downweight the signal 
close to the cluster core where cluster dwarf galaxies may pollute 
the signal and where the reduced shear g rather than the shear y 
should be used. The parameter in Eq. (O controls the shape 
of the filter function. Eor low values of x^, more weight is put to 
smaller filter scales. In the absence of the exponential cut-offs, 
variations in x^ and flmax are obviously degenerate, but the ex- 
ponent ial cut-off makes the weight function (|9]l a two-parameter 
family. iHetterscheidt et al.l (j2005) find that the optimal value of 
Xc = 0.15 for a range of cluster masses and redshifts, so we will 
use this value throughout this work. 

The projection of large-scale structure (LSS) along the line 
of sight can potentially be a serious contaminant for every weak 
lensing observation of galaxy clusters. Such projections of sheets 
and filaments inevitably lead to false cluster detections at all sig- 
nificance levels expected from real clusters dHSOSh an d mis sed 
cluster detecti on exce pt for the most massive systems (|H04|) . It 
was shown bv lHS05l that, even in the absence of shape noise, 
the efficiency of a weak lensing search for galaxy clusters does 
not exceed ~ 85% and also depends on the shape of the filter 
fun ction. 

iMaturi et alj (l2005b propose to treat the projections of 
masses in the background and foreground of the cluster, i.e., the 
large-scale structure as a source of non-white noise. They con- 
structe an optimal filter that tries to maximize the signal caused 
by a galaxy cluster while downweighting the cosmic shear sig- 
nal on scales of interest. The resulting filter function Qlss de- 
pends not only on the expected shear profile of the galaxy cluster 
but also on the number density of EBG and on the convergence 
power spectrum to describe the expected noise. The latter noise 
contribution in turn depends on the redshift distribution of the 
EBG. In the case of pure white noise, the Maturi filter takes the 
shape of an NEW shear profile, i.e., essentially the same form 
as the Schirmer filter, and both filters are virtually equivalent. If 
projections of the LSS are a significant noise source, then using 
2lss instead of Qnfw can lead to significantly lower contamina- 
tions of the cluster catalog with false positives. 

5. 1 . Signal and noise of fhe aperture mass 

Before we apply the aperture mass statistic to the XES data we 
need to understand the properties of the aperture mass function 
in some more detail. Equations (|5]l and (|6]l allow us to compute 
the SNR of Map as 



(10) 



where Qi is the weight assigned to the ith galaxies by the radi- 
ally symmetric weight function This expression, however, 
considers only the noise caused by the random ellipticities and 
not the shot noise of the finite sampling of the Map statistic by 
the population of background galaxies. 

We used ray-tracing simulations to examine these effects 
in detail. These were generated by tracing light rays through 
i Jie ACDM y-bod y simulation of the VIRGO consortium 
jjenkins et alJll998l) . These A^-body simulations have been car- 
ried out with the following parameters: Qm - 0.3, Qa - 0.7, 
h-iQ - 1, (Tg = 0.9, E = 0.21, and the index of the primordial 
power spectrum Us - I. The population of background galaxies 
has been assumed to be a 5-function peaked at z = 1 . 
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We now give a short description of our ray - tracin g algo- 
rithm. More details of can be found in 'Hartlagl ( l2005h . From 
redshift to the source redshift, 1024^ light rays are traced 
through 16 slices of 202.9 /i^,} Mpc thickness onto an output grid 
of 1 X 1 square degree. Each redshift slice corresponds to one 
output box of the A^-body simulation and is projected as a whole 
onto a lens plane, preserving the periodic boundary conditions of 
the A^-body box. To avoid repetition of structure along the line 
of sight, the planes are randomly shifted and rotated. The light 
rays are shot from the observer through the set of lens planes, 
forming a regular grid on the first plane. We then use FFT meth- 
ods to compute the lensing potential on each lens plane, from 
which we obtain the deflection angle and its partial derivatives 
on a grid. The ray position and the Jacobian of the lens mapping 
for each ray are obtained by recursion. Given the ray position on 
the current lens plane, its propagation direction (known from the 
position on the last plane), and the deflection by the current plane 
interpolated onto the ray, we immediately obtain the ray position 
on the next plane. Differentiation of this recursion formula with 
respect to the image plane coordinates yields a similar relation 
for the Jacobian of the lens mapping, taking the tidal deflection 
field computed before into account. The recursion is performed 
until we reach the source plane. From the final Jacobian, we ob- 
tain noiseless convergence and shear maps in the usual way. 

Since we wish to study which convergence/Map peaks cor- 
respond to real dark matter halos, we also have to compute the 
lensed positions of the central particles of the dark matter halos 
contained in the halo catalogs. We achieve this by projecting the 
halo position onto the lens planes and identifying the light ray 
that passes closest to the halo. The lensed position of the halo is 
then computed by inverting the linearized lens mapping around 
this ray. 

Fifty diff'erent realizations were made by using in each case 
different random shifts and rotations of the lens planes. Lensing 
catalogs were created by randomly distributing galaxies with 
an ellipticity dispersion of cr^ = 0.38 over the output grid of 
the ray-tracing simulations on areas corresponding to the sizes 
of actual XFS fields until the number density of the respective 
XFS field was reached. When placing galaxies, the masks used 
in the catalog generation of the real data were applied to simu- 
lated catalogs as well, to realistically model the effect of holes 
in the field. Three different masks and number densities were 
taken from the XFS data. All masks were put on each of the 50 
ray-tracing realization to obtain a high number of lensing sim- 
ulations. The effective (unmasked) area and number densities 
of these fields are 970arcmin^ and 13.4arcmin"^, 1022arcmin^ 
and 17.7arcmin"^, and 981 arcmin^ and 18.2 arcmin"^. The total 
area covered by our ray-tracing simulations is 41 .3 sq. deg. 

With typical sizes of the XFS fields of 35' x 35', the indi- 
vidual lensing simulations are not totally independent but have 
some overlap because the side length of one ray-tracing simu- 
lation is only 60', not enough to accommodate two XFS masks 
next to each other. The catalogs were placed on the ray-tracing 
simulations such that this overlap is minimized. While not com- 
pletely independent, these overlapping areas were covered by 
catalogs with different maskings, different number densities, and 
different realizations of Gaussian noise so that for our purpose 
- understanding the noise properties of Map from realistic sim- 
ulations - no significant correlation between individual lensing 
simulations is expected. 

We computed the aperture mass for 9 different filter scales 
and the filter funct ions proposed by Schirmeret al. (2007) and 
iMaturi et al.l (l2005l) from the same catalogs. As both filter func- 
tions are based on an NFW model of a cluster, we chose the filter 



Table 2. Filter radii for Map computation and corresponding 
virial mass. 



Radius 


Mvi,- 


(kpc//i7o) 


(10'^/j,-;Mo) 


1000 


0.76 


1247 


1.5 


1493 


2.5 


1740 


4.0 


1986 


6.0 


2233 


8.5 


2479 


12.0 


2726 


15.0 


2972 


20.0 



radii based on an assumed fiducial cluster model. In this model 
the cluster is at redshift z - 0.3, the redshift at which we expect 
the lens strength in our survey to be maximized. Table |2] gives 
these filter radii and the corresponding virial mass if the filter 
scale is interpreted as the virial radius of the cluster We also 
need to fix the redshift of the source galaxies to model the large- 
scale structure power spectrum in the Maturi filter We assumed 
that all background galaxies are at Zs - 0.8 and that the number 
density of background galaxies is n = 18arcmin"^. We com- 
puted Map on a grid and set the pixel size in this grid such that 
one pixel corresponds to 50 h^^ kpc at the redshift of the fiducial 
cluster model. 

For all lensing simulations and filter scales we computed Map 
and -Map (for later detection of negative peaks as a control) for 
both filter functions from the input catalog and Map after rotating 
(1) all galaxies by 45°; (2) every galaxy by a random angle for 
the Schirmer filter. Maps of the aperture mass and their SNR 
(S-maps), which were computed as well, were saved as FITS 
images. 

As for the detection of peaks in the matched filter maps, we 
used SExtractor to identify shear-selected peaks in the Map- 
maps. For this purpose we run SExtractor in dual-image mode 
with the Map-map as the detection image and the S -map as the 
measurement image. This means that the SNR of a Map peak 
is computed with Eq. ( fTOt and not determined by SExtractor. 
This is a small differenc e fro m the matched filter pipeline, in 
which we did not use the |P96| likelihood to determine a signifi- 
cance but used SExtractor detection significances. 

The detection threshold is set to 2cr and the minimum de- 
tection area scales with the filter such that it corresponds to the 
pixels covered by a circle with a radius twice as large as the 
filter scale. SExtractor is run without deblending, i.e., every 
contiguous area above the detection threshold is counted as one 
object. Peaks found in this way in different filter radii were as- 
sociated based on positional coincidence. 

We first examine the results obtained with the Schirmer fil- 
ter before comparing these to the ones obtained with the Maturi 
filter Positional offsets between the weak lensing peak positions 
and the position of the BCG or the center of X-ray emission from 
the ho t intracluster gas are commonly o bserved (IWittman et al.l 
'2006*) and expected Dietrich et all2005h . The siz e of flie offsets 
has been studied for the case of an isolated SIS bv lDietrich et al.l 
(i2005b . The ray-tracing simulations used here allow us to in- 
vestigate the additional effect of large-scale structure along the 
line-of-sight. 

The catalogs of weak lensing halos produced in this way 
from the Map images were associated with a catalog of actual 
dark matter halos in the VIRGO simulation. For this purpose 
we only considered halos with masses in excess of lO'^^M© and 
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Fig. 3. Cumulative distance distribution of the 434 Map peaks 
that could be associated to a dark matter halo in the VIRGO 
simulation. 75% of all matches are made within a 2;i5 radius. 



with redshifts 0.1 < z < 0.7 as these are roughly the ones to 
which we expect to be sensitive in our galaxy cluster survey. 
The maximum distance allowed for a match between halo po- 
sition and Map peak was determined by the virial radius of the 
halo and the size of the Map peak as determined by SExtractor. 
Note that the position of the halo was derived from the most- 
bound particle in that halo identified by the friend-of-friend halo 
finder employed to generate the halo catalog. In rare cases this 
is in the center between what by eye be would identified as two 
separate halos. 

In the ray-tracing simulations, 434 peaks in the Map-maps 
could be associated with dark matter halos in the VIRGO sim- 
ulation. Figure |3] shows the cumulative distribution of their po- 
sitional differences. From the number density of Map peaks and 
the average size of the association radius we estimate that ~ 100 
or ~ 25% of those matches are chance coincidences. Inspection 
of Fig.[3]shows that 75% of all positional offsets are smaller than 
2! 15, which is the maximum offset we will allow from here on. 
Note that on the one hand this is significantly smaller than the 
3' matching radius adopted by HS05, who used a higher number 
density and a smaller ellipticity dispersion. On the other hand, it 
is significantly larger than the offsets found for an isolated SIS by 
iDietrich et al.l (l2005l) . This possibly indicates a non-negligible 
influence of large-scale structure along the line of sight on the 
weak lensing peak positions of dark matter halos. One, however, 
has to be careful when drawing this conclu sion as we are looking 
at halos at very different redshifts, while iDietrich et al.1 (|2005[) 
studied only systems at one redshift. 

Aperture mass peaks not associated with dark matter halos 
can be caused either by projections of large-scale structure mim- 
icking a shear signal of a cluster or by the shape noise of back- 
ground galaxies that can cause random tangential alignments. In 
real data the measured ellipticities must be corrected for atmo- 
spheric seeing and PSF distortions of the instrument. Residuals 
in this correction can lead to spurious alignments of background 
galaxies, both in curl-free shear fields (pure E-modes) and in 



non-curl free shear fields (also including B-modes). E-modes are 
transformed into B-modes (and vice versa) by rotating all galax- 
ies by 45°. Because gravitational lensing only creates E-modes, 
the observed power of B-modes is often used as a quality check 
of the PSF correction. We did not check our lensing signal for the 
presence of B-modes. However, cosmic shear studies using PSF 
correction schemes very similar to ours find that B-modes typi- 
cally occu r on scales smaller than are relevant for our Mao ker - 
nels (e.g.. IVan Waerbeke et al]|2005l: iHetterscheidt et al.ll2006l) . 
It is thus safe to assume that we can apply the results of our 
ray- tracing simulations to the XFS. 

Figure|4]shows the significance distribution of shear-selected 
peaks in the different kinds of aperture mass maps created with 
the Schirmer filter These are (1) the distribution of peak signif- 
icances for all Map peaks created from the ray-traced catalog, 
i.e., those peaks one would find in real data, (2) negative peaks 
or, considering how the peak finding pipeline is run, peaks found 
in -Map-maps, (3) peaks found in B-mode Map-maps, i.e., maps 
of aperture mass created after all galaxies in the ray-traced cat- 
alog were rotated by 45°, (4) weak lensing peaks that could be 
associated with dark matter halos within a matching radius of 
2^15, and (5) peaks found in mock catalogs with random elliptic- 
ities. It is important to emphasize that the random seed was kept 
fixed so that when computing Map on different filter scales the 
orientation of galaxies in the input catalog remained unchanged. 
Here we consider only those peaks that were detected on at least 
two filter scales. We justify this choice later when examining the 
influence of the number of filter scales in which peaks are de- 
tected. The significance used in the histogram is the maximal 
significance values of all filter scales a peak was detected in. 

A number of features in Fig. |4] are worth a more detailed 
discussion. First, we note that the number density of peaks in 
the Map- and -Map-maps declines towards lower significances, 
while in the B-mode and random maps it remains roughly con- 
stant below 3.25cr. Naively, one would expect an increase in all 
curves towards lower significances. The observed behavior is 
due to a selection bias when running SExtractor on the FITS 
images. The detection threshold is derived from the standard de- 
viation of the background in these images. The B-mode and ran- 
dom maps are overall flatter than the E-mode Map-maps. This 
leads to a detection threshold at which peaks of lower signifi- 
cance in the S-maps are detected in the Map-maps for random 
maps and B-mode maps than for E-mode maps. 

Second, the number densities of random and B-mode peaks 
are very similar, with the latter slightly lower As we did not sim- 
ulate the systematic influence of instrumental PSF corrections, 
the only sources of B-modes are the shape noise of background 
galaxies and finite field effects in the Map estimator. The effect of 
shape noise alone is simulated by the random catalogs, while the 
B-mode peaks are a combination of shape noise and systematic 
effects due to the finiteness of the field and holes in the data. The 
fact that the number density of B-mode peaks is compatible with 
the number density of random peaks shows that systematic ef- 
fects contributing to B-modes are not an important noise source. 
This is confirmed by a visual inspection of the peak distribution 
indicating that B-mode peaks do not show an obvious tendency 
to appear close to holes in or edges of the data field. 

Third, the number of E-mode peaks is higher than the num- 
ber of any other peak statistic in all significance bins > 3.25cr 
in which the aforementioned selection plays no role. The sum 
of the true halos peaks and the random peaks is compatible with 
the number of E-mode peaks in the significance bins from 3.25cr 
to 4.25cr, if one assumes Poissonian statistics. At higher signifi- 
cances, an excess of E-mode peaks is observed. We surmise that 
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Fig. 4. Peak significances of shear selected peaks in the sim- 
ulated maps of Map (black solid), -Map (green long dashed), 
B-mode Map (red short dashed), and mock catalogs (blue dot- 
dashed). The solid pink line corresponds to peaks successfully 
associated with dark matter halos within 2'.15. The thin light- 
blue dotted line is the sum of the random peaks (blue) and the 
true halos (pink). 



Fig. 5. Number density of lensing peaks in dependence on the 
number of filter scales in which the peak is detected. The 
color/line-style coding is the same as in Fig. |4] Dark matter 
halo peaks show virtually no dependence on the number of fil- 
ter scales, while the number density of spurious peaks sharply 
declines with a requirement on the minimum number of filter 
scales. 



this is due to projections of large-scale structures. This means 
that we can expect a significant fraction of spurious peaks at 
almost all s ignifi cances, a result that is compatible with earlier 
findings of |ho3 and|HS05. At low significances the spurious 
peaks will be dominated by shape noise, while at high signif- 
icances many spurious peaks will be caused by the projection 
of large-scale structures. We emphasize that the latter class of 
peaks is in fact caused by gravitational lensing. They just do not 
correspond to a single mass concentration in 3-d space. These 
peaks are spurious peaks only in the sense of a galaxy cluster 
search. 

Fourth, negative peaks are relatively rare. The two effects 
leading to this result are best understood in terms of the filter 
function U that is related to Schirmer's 2nfw function by Eq. (|4]i 
and acts on the surface mass density. The function U has a nar- 
row positive peak close to the origin with extended and shal- 
low negative wings to satisfy the condition that U has zero total 
weight. The comparably low number density of -Map peaks is 
then caused by (1) the shallowness of negative wings, which will 
limit the peak strengths of negative peaks; and (2) by U acting 
as a bandpass filter for structures with the same size as the char- 
acteristic filter scale. The large extent of the negative wings will 
make negative peaks more extended than positive peaks, natu- 
rally leaving less space for other peaks. 

By associating peaks found in different filter scales with each 
other, we can also examine whether the number of filter scales 
Hf that a peak is detected in says something about the correspon- 
dence of the lensing signal to a dark matter halo. Figure|5]shows 
the number density of peaks detected in the ray-tracing simula- 
tions as a function of rif. The color/line-style coding is the same 
as in Fig.m One clearly sees that the number density of lensing 



peaks associated with a dark matter halo is virtually independent 
of the number of filter scales the Map peak appears in. Real dark 
matter halos show up as often in only 3 filter scales as they do in 
9 filter scales. 

The behavior of spurious peaks is very different. Their num- 
ber falls off monotonically as a function of «f , with an exception 
in the last bin, i.e., the peaks that are detected in all filter scales. 
This dependence can be used to impose a selection criterion to 
decrease the contamination of shear selected clusters with spuri- 
ous lensing peaks by requiring that a Map peak must occur in at 
least a given number of filter scales. However, as stated above, 
this will always exclude a number of real clusters as well. As a 
compromise between efficiency and completeness, we imposed 
the condition that «f > 3 in the XFS. 

We now briefly compare these results to the ones obtained 
with the Maturi filter on the same catalogs. This filter can poten- 
tially result in fewer spurious peaks and increase the SNR of real 
clusters, thereby increasing the number of real cluster detections. 
We find that the results of the Schirmer and Maturi filters are en- 
tirely consistent for our survey parameters. The number of Map 
peaks associated to dark matter halos within the 2! 15 matching 
radius (316) is, contrary to expectations, slightly smaller than 
for the Schirmer filter (325). The difference might well be due 
to a lower contamination with spurious peaks at high signifi- 
cances where the contamination with LSS projections should be 
suppressed by the Maturi filter. The observed peak offsets from 
the dark matter halo positions are consistent with the those of 
the Schirmer filter. We compared the peak significances of Map 
peaks related to dark matter halos computed with both filters and 
find {S(Qlss) -S(Qt^Fw)} = -0.06 ± 0.27, i.e., the peak signifi- 
cances are consistent with each other, non-significantly favoring 
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the Schirmer filter. Again, this small difference might be due to 
the relative suppression of highly significant projections of the 
LSS. These findings agree with Maturi et al. (2007) who com- 
pare the performance of their filter with the cluster search by 
ISchirmer et al.l (|2007|) on the same fields. As a result of our com- 
parison we limit the weak lensing cluster search in the XFS to 
the Schirmer filter that is less complex and faster to compute. 

5.2. Weak lensing catalog 

Based on the optical catalogs created with SExtractor's default 
convolution kernel (Sect. \4.1\. we created lensing catalogs. We 
used the KSB algorithm CKaiser et al.lll995 ) to obtai n shear esti- 
mates closely following the procedure described by lErben et al.l 
(2001). From the KSB catalogs we constructed catalogs of prob- 
able background galaxies that are reliable shear estimators by 
imposing the following selection criteria: Objects with SNR< 2, 
Gaussian radius < 0'.'33 or > I'.' 19, or corrected ellipticity 
e > 1.0 were deleted from the catalog. We also deleted objects 
whose pre-seeing shear polarizability tensor P'' has a negative 
trace, and bright galaxies with R < 21 mag. 

Table[T]gives the effective (unmasked) area of all XFS fields 
used for the weak lensing cluster search, as well as their num- 
ber density of galaxies in the weak lensing catalogs. The total 
area used for weak lensing is 6.4 sq. deg. The number density of 
background galaxies averaged over this area is 14.1 arcmin"^. 

The aperture mass statistic with the Schirmer filter function 
was estimated from the resulting catalogs on the filter scales 
listed in Table |2l No weighting of individual galaxies was done. 
Peaks in the M^p-maps were detected as described in the previ- 
ous section. The final catalog lists all aperture mass peaks with 
maximum SNR > 5, or with maximum SNR > 3 if the peak 
has an X-ray or matched-filter counterpart, or was previously re- 
ported as cluster (candidate) in the literature within 2^15 radius 
from the lensing position. All shear peaks must be present in at 
least 3 filter scales. 

In the following we give comments on some specific fields. 

- T Leo - No association could be found for the 4 Map peaks 
with CTtnax > 3 and «f > 2 in this field. Thus, none of these 
peaks are included in our catalog. 

- Field 864-1 - None of the 5 Map peaks with cr^ax > 3 and 
Hf > 2 could be associated with X-ray or matched-filter clus- 
ters. These peaks are not included in our catalog. 

- Field 864-9 - This field contains one of the three weak- 
lensing selected cluster candidates that was previously iden- 
tified in the literature as a cluster candidate but not found 
by the matched filter or in X-ray. The candidate BLOX 
J1343.5-0022.8 is outside the FOV of XMM-Newton. The 
cluster candidate is very elongated and might for this reason 
be missed by the matched filter pipeline. 

- LBQS 2212-1759 - This field contains the remaining two 
weak-lensing cluster candidates that match only clusters pre- 
viously reported in the literature, but were found with neither 
the matched filter nor the X-ray survey. Both cluster candi- 
dates match optically selected clusters from '007*, see also 
Sect. 14.21 Gavazzi & Soucail (2007) have performed a weak 
lensing cluster search on this field as part of the CFHTLS 
Deep fields and have not found any convergence peak in the 
FOV of WFl above their detection threshold of 3.5cr. 

We selected a total of 3 1 cluster candidates using the aper- 
ture mass method on 23 WFl fields. The full catalog of shear 
selected cluster candidates is available in electronic format at 



the CDS. There we list the position of the Map peaks, the signifi- 
cance of their detection, the number of filter scales a cluster was 
detected in, and the filter scale in which the SNR of a peak was 
maximized. This catalog contains 4 of the 5 clusters that were 
primary XMM-Newton targets in the XFS fields. The fifth clus- 
ter, MKW 9, is at redshift z = 0.04, too low to be detectable with 
weak lensing. 

On average we find 1.3 weak lensing cluster candidates per 
XFS field. Most Map peaks do not correspond to a cluster can- 
didate. This is to be expected from the results of our ray- tracing 
simulations and illustrated by Figs.|4]and |5] The two fields con- 
taining only unmatched Map peaks are statistically expected. The 
absence of weak lensing peaks in these fields does not hint at 
problems with their lensing catalogs. 

6. Summary and discussion 

In this work we have selected galaxy cluster candidates indepen- 
dently with three different methods: optical matched filter algo- 
rithm, extended X-ray emission, and the shear signal induced by 
massive foreground structures. 

We found a total of 155 cluster candidates in 23 WFl fields, 
or 24.2 cluster candidates per square degree. Most cluster candi- 
dates were found with the optical matched filter (116), f ollowed 
by X-ray emission (59). As was previously shown (H04; HS05|) 
and confirmed by our lensing simulations using ray-tracing sim- 
ulations (Sect. 15. Il l, the efficiency of weak lensing for cluster 
selection is relatively low. To avoid being dominated by spuri- 
ous weak lensing signals, we limited the catalog of weak lensing 
selected cluster candidates to those that have an optical or X-ray 
counterpart, either found in our own survey or previously re- 
ported in the literature. We found significant lensing signals for 
31 cluster candidates, of which 12 are previously known cluster 
(candidates). Eleven of the weak lensing selected clusters were 
detected with both the matched filter and X-ray emission, ex- 
cluding A 1882 that is not part of our X-ray catalog; 6 of these 
are previously unknown cluster candidates. 

Comparing the redshift estimates of the Postman matched- 
filter method to spectroscopically measured redshifts or other 
photometric estimates, we find that these work surprisingly well. 
The mean difference in redshifts (Az/(1 + Ztme)) is marginally 
significant only for the spectroscopic sample with -0.06 ± 0.05, 
and consistent with zero if we also trust the redshifts of the 
photometric sample that gives (Az/(1 + Ztme)) - 0.06 ± 0.07. 
Considering that we use only one passband to derive the redshift 
of matched-filter clusters, this result compare s favorably to wha t 
is achievable with more colors. For example. iGoto et all (l2002h 
report a mean Az for their "cut-and-enhance" method of 0.02 us- 
ing four colors, but only after outliers with Az > 0. 1 have been 
rejected. However, the matched filter estimated redshifts come 
with a significant bias, which puts higher redshift cluster at red- 
shifts that are typically too low by Az = 0.1 . . . 0.2. 

We described in detail how we developed the selection cri- 
teria for our weak lensing sample using ray-tracing simulations. 
We find that - at least for the comparatively low number densi- 
ties and peak significances we are dealing with - the dominant 
source of noise is the shape noise of the background popula- 
tion and not projections of the large-scale structure. Our ray- 
tracing simulations suggest that the contamination with projec- 
tions of the large-scale structure becomes more important at 
higher significances. However, the area covered by the XFS is 
comparatively small and the absolute number of highly signif- 
icant Ma p peaks is consequ ently small. The Map kernel devel- 
oped bv iMaturi et al.l (l2005l) to minimize the effect of large-scale 
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structure on weak lensing cluster searched thus could not per- 
form b etter than the filter function proposed by ISchirmer et alj 
(I2007h used in this work. It is sensible to assume that the Maturi 
filter will perform better on deeper surveys dPace et al.ll2007h for 
two reasons: (1) The higher number density of FBG will reduce 
the white noise component and (2) a deeper survey will probe a 
larger volume and hence increase the contribution of LSS pro- 
jections. 

The matching radius of 2! 15 we deduced from our ray- 
tracing simulations is considerably smaller than the 3' employed 
by[HS05. Considering that HS05 used a higher number density 
of FBG with an ellipticity dispersion lower by a factor of 1 .6, 
even their already low efficiency estimates for weak lensing clus- 
ter searches still seem to be too optimistic. 

The number density of weak-lensing selected cluster candi- 
dates is 4.8 per square degree. This is slightly lower than the 
number density of Map peaks with a halo counterpart in the 
ray-tracing simulations, which is 6.1/sq. deg. However, the av- 
erage number density of background galaxies in the XFS data 
is only 14.1 arcmin"^; this is somewhat lower than in the ray- 
tracing simulations, which had an average number density of 
16.5 arcmin"^. Whether the difference in cluster counts can re- 
ally be attributed to the difference in number density should be 
checked by adjusting the simulation parameters to match the 
XFS observations better. The trend to slightly lower number den- 
sities is also present if we select only Map peaks with a higher 
SNR> 4. The XFS contains 1 1 of these highly significant peaks 
associated with a matched filter or X-ray cluster candidate. This 
corresponds to 1.7/sq. deg., compared to 2.3/sq. deg. in the ray- 
tracing simulations. These high significance peaks include all 
4 clusters that were primary XMM-Newton targets at redshifts 
accessible by weak lensing. These 4 clusters alone contribute a 
number density of 0.6/sq. deg. 

With the cluster sample presented here we have built a 
solid foundation for studying possible selection effects in either 
method. This will be the subject of a follow-up paper (Dietrich 
et al., in preparation). 
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Appendix A: X-ray observations 

Table I A. O gives a summary of the EPIC X-ray observations contributing to the 
XFS data used in this work. The table gives for each field: in Col. 1 the field 
name; in Col. 2 the XMM-Newton observation ID; in Col. 3 the nominal ex- 
posure time in seconds; in Col 4-6 the settings for each of the cameras. Here 
(E)FF indicates (extended) full frame readout, LW large- window mode, and SW 
small-window mode. These cameras and their settings are described in detail in 
Ehle et al. (2004). For some fields additional observations were available, but 
these were discarded mainly due to unsuitable camera settings. 

"Field 864" is a mosaic of 3 X 3 XMM-Newton observations. The name was 
assigned by the principal investigator Fields 864-2, 4, 5, 6 are listed in TablelAJ] 
although they are not covered by our WEI observations. The reason is that these 
data overlap with the Fields 864-6, and 9, which are XFS survey fields. 
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Table A.l. Information about X-ray imaging in the XFS. 



Field 


Obs. ID 


^exp / S 




Camera settings 




BPM 16274 


0125320401 


33 728 


EPN EFF 


MOSl FF 


M0S2 FF 




0125320501 


7845 


EPNFF 


MOSl FF 


M0S2 FF 




0133120301 


12 022 


EPN FF 


MOS 1 FF 


M0S2 FF 




0133120401 


13 707 


EPN FF 


MOSl FF 


A iff /~V O 1 <| ' 

MOS2 FF 




01/5320/01 


45 951 


brJN rr 


MUol rr 


AT/~\co trtr 
MUSz rr 




01 j3yj0101 


jljD 


COM T7T7 


MUo 1 rr 


A 4/^ CO T7T7 

MUSz rr 




UU41 1 /UlUl 


^ 1 TO/1 

J 1 /Z4 


nrJN r,rr 


JVlUo 1 rr 


MUiSz rr 


RX J0505. 3-2849 


0111160201 
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EPN EFF 


MOSl FF 


MOS 2 FF 


RBS 0864 


0108670101 


56459 


EPN FF 


MOS 1 FF 


M0S2 FF 


QSO BO 130-403 


01 12630201 


37 870 


EPN FF 


MOSl FF 


A It /~\ O ^ 1 <| < 

MOS2 FF 


D r) 1 1 Aim 


U1j0o/U4U1 


31 418 


bPN FF 


MUSI rF 


MUSz FF 


SDSS J 104433.04-012502.2 


0125300101 


62 310 


EPN FF 


MOSl FF 


A iff i~\ O 1 '1 ' 

MOS2 FF 


MS1054.4-0321 


AAA A OAA 1 A t 

0094800101 


41 021 


EPN FF 


MOSl FF 


A iff /~V O 1 '1 > 

MOS2 FF 


T Ti"" 1 1 n /I 1 one 

HE 1104-1805 


Al I'^/'TAIAI 

01 12630101 


36428 


EPN FF 
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A iff /~\ O ^ 1 <| ' 

MOS2 FF 


PO 1115+080 


AAO'^ O /I A 1 A 1 
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MOS 2 FF 
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iVlUoZ rr 


iNOC 40 DO 


Ul iUVoUZUi 


CO OTO 

JO Zj / 


nriN nrr 


IVlvJo i rr 


iVlL/oZ rr 


yjpyj D iZ4D— Uj / 


UUuUj /UzUl 


A 1 OTQ 

41 11 


crJN rr 


JVlUo 1 rr 


A /toco cc 
MUoz rr 


T^Ial^l OA/1 1 

rielcl oo4-i 


A 1 1 1 OO 1 AA 1 


1 r\ "ITT 
10 J / / 


nrJN nrr 


MUSI rr 


MUSz rr 


Field 864-2 


0111282401 


IQll 


EPN EFF 


MOSl FF 


A Iff O ^ 1 1 1 « 

MOS2 FF 


Field 864-4 


0111281301 


14541 


EPN EFF 


MOSl FF 


M0S2 FF 


Field 864-5 


0111281401 


8643 


EPN EFF 
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EPN EFF 
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M0S2 FF 
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8623 


EPN EFF 


MOS 1 FF 
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23 567 
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MOS2 FF 


MKW 9 


0091140401 


45 414 


EPN EFF 


MOSl FF 


A Iff O 1 * 1 < 

M0S2 FF 


LBQS 2212-1759 


0106660101 


60508 


EPNFF 


MOSl FF 


M0S2 FF 




0106660201 


53 769 


EPNFF 


MOSl FF 


M0S2 FF 




0106660401 


35114 




MOSl FF 


M0S2 FF 




0106660501 


11568 


EPNFF 


MOSl FF 


MOS2 FF 




0106660601 


110 168 


EPNFF 


MOSl FF 


MOS2 FF 


NGC 7252 


0049340201 


28 359 


EPNFF 


MOSl FF 


M0S2 FF 


PHL 5200 


0100440101 


46 681 
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MOSl FF 


MOS2 FF 
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